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Abstract. EE Cep is an eclipsing binary with a period of 5.6 years. The next eclipse will occur soon, in May-June 
2003, and all available past eclipses were collected and briefly analysed. EE Cep shows very large changes of 
the shape and the depth of minima during different eclipses, however it is possible to single out some persistent 
features. The analysis suggests that the eclipsing body should be a long object surrounded by an extended semi- 
transparent envelope. As an explanation, a model of a precessing optically thick disc, inclined to the plane of 
the binary orbit, is invoked. The changes of its spatial orientation, which is defined by the inclination of the disc 
and the tilt, induced most probably by precession of the disc spin axis with a period of about 50 years, produce 
strange photometric behaviour of this star. The Ha emission, and possibly the Nal absorptions, show significant 
changes during several months outside of the eclipse phase. 
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1. Introduction 

EE Cep (BD+55°2693) is an extraordinary long-period 
eclipsing binary showing very large and complex changes 
of the shape and the depth of the minima (e.g. Baldinelli 
& Ghedini 1977). Quite recently Mikolajewski & Graczyk 
(1999]- hereafter MG) have suggested that the reason for 
these changes is the presence of a non-stellar eclipsing 
body in the system - a precessing, opaque and dark disc. 

This star was identified as a variable star by Romano 
( 1956) in 1952 and it was confirmed by Weber (195HJ) who 
reported observations of the previous minimum in 1947. 
The colour changes during the 1952 minimum were very 
small and the star was originally classified by Romano 
as an R CrB variable. Perova (1957) reported about 500 
photographic observations of EE Cep between 1899-1908 
and 1933-1955 but no brightness changes were recorded. 
Romano & Perissinotto £1966) have suggested that it may 
be an eclipsing binary but they did not give any approx- 
imate ephemeris. Its binary nature was confirmed later 
by Meinunger l|1973fl who analysed over 500 photographic 
plates from Sonneberg Observatory. The period of the star 
is almost 2050 days and only the primary minimum can be 
detected. Meinunger has not detected any changes of the 
brightness between the 1958 and 1964 minima. The bright- 
ness of the star outside of the eclipse phases is constant 
within the observational errors: V max — 10.78 (Meinunger 
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& Pfau ll9cH1 MG). The primary component is apparently 
a normal B5 III star (Herbig 113501 Baldinelli et al. lTTMjl . 
The emission visible in H a and Hp lines arises most likely 
from an invisible companion (Bruckner 1976 ). The absorp- 
tion spectrum of the companion is never seen. 

The early model of EE Cep presented by Meinunger 
(1976) assumed that the companion should be a cool star 
of M type, much less luminous but larger than the primary 
B5 III star, surrounded by a strongly extended envelope. 
The pulsations of the cool star would change its diameter 
and thus lead to the changes of the minima duration and 
depth. The envelope would cause an atmospheric eclipse in 
the beginning and in the end of the minimum, producing 
typical wings in the light curves. However, such a simple 
model cannot explain some of the features observed in EE 
Cep, as was argued by MG and in this paper. 

2. Photometric behaviour 

2.1. A review of past eclipses 

The typical minimum of EE Cep has large wings, the du- 
ration of the eclipse is D sa 40 days or a phase inter- 
val AP = 0.02, and the depth of the eclipse is moderate 
AB Bs 1.5 mag. Such minima were observed in 1975, 1980 
and 1997. But some eclipses observed up to date devi- 
ated strongly from this picture. We have collected obser- 
vations of all previously observed eclipses of EE Cep in 
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Fig. 1. EE Cep's eclipses in B/m pg light. 



two roughly similar photometric wavelength bands: 1) B 
and m pg ; 2) V, m pv and m v - ls . The references used for the 
photometry are listed in TablcQ whereas the data are pre- 
sented in Figs. ^ an d [21 Part of the historical data were 
published only as diagrams and they needed to be read 
from these figures. All observations were phased using the 
ephemeris given by MG. 

The 1947 minimum {E = —1) was the first one de- 
tected and its depth was 1™5 in m pg . There is no infor- 



mation about the duration and shape of the eclipse and it 
is not shown in Fig. ^ The 1952 minimum (E = 0) was 
deep (at least 1.7 mag in m pg ) and has an asymmetric 
shape. During the next eclipse, in 1958 (E = 1), there 
are only two observational points and two upper limits 
« 13"Mn rripg at the bottom of the minimum (Romano 
& Perissinotto 1966), therefore we also omitted this min- 
imum in Fig.^ Nevertheless, it could be very similar, or 
maybe a little deeper, to the following one in 1964 (E = 2). 
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Fig. 2. EE Cep's eclipses in V/m pv /m vis light. 



The depth of these two minima is at least 2 mag. The du- 
ration of the 1964 eclipse was only 30 days. 

During the 1969 (E = 3) eclipse, the minimum was, in 
contrast to the previous one (E = 2), extremely shallow 
(Am pg = 0.6 mag) and a flat bottom phase was clearly ob- 
served. Also, the 1969 eclipse was probably the longest one 
- its total duration was about 60 days. The 1975 (E = 4) 
and 1980 (E = 5) eclipses were more or less typical, while 
the 1986 {E — 6) minimum was again shallow and strongly 
asymmetric. The 1992 (E = 7) eclipse was observed visu- 



ally by three amateur astronomers from AAVSO and seven 
VR measurements with a CCD camera have been ob- 
tained by J. Borovicka from Ondrejov Observatory in the 
Czech Republic (Halbach 1992). Unfortunately, there is a 
zero-point offset of +0™16 between Borovicka's V magni- 
tude scale and Meinungers's l|1976f) and MG's scales. Also, 
the visual light curves strongly suffer from differences in 
zero-point scales. Thus we have shifted visual estimations 
to the corrected CCD observations, although three points 
still lie out of the light curve (Fig. [21 . This minimum was 
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Table 1. The duration and the depth of eclipses. 



Epoch 


Year 


AP 


A V/?77,.rw. / TYl,n. n c 
i — i V 1 1 1 vpv j I I t"uiS 


A R /rn ir .r, 

. — \ i y j 1 1 vpn 


References (source of data,) 




1 Q47 






L.O. 


VVcUcl ^juU| ^ LcLiJlc ) 


n 

u 


1 

J- lyO£ 






1 7+n 9 




I 


1958 


01 5+0 003 




1.9±0.2 


P? mn p n n At" Ppvi qqi n f n 1 Qf"i(~i ftnnlpi 


2 


1964 


0.015±0.002 




2.1±0.2 


Meiirunser jl£)73 f fi&'ure) 


3 


1969 


0.030±0.003 




0.6±0.1 


Baldinelli et al. H17751 (table) 


4 


1975 


0.021±0.001 


1.34±0.02 


1.43±0.02 


Meinuneer |I976| ( table) . Zaitseva et al.|I975| ("table) 


5 


1980 


0.019±0.001 


1.25±0.02 


1.30±0.02 


Baldinelli et al.flMIlffieure). Meinuneer & Pfau[T9"8"Tl (table) 


6 


1986 


0.025±0.002 


1.1+0.1 




Di LucaEnHH (figure) 


7 


1992 


0.011±0.002 


0.85±0.15 




Halbach IT^l (figure) 


8 


1997 


0.023±0.001 


1.61±0.02 


1.68±0.02 


MG, Halbach 1999 (figure) 



Note: the total duration was calculated as time period between l"st and 4+th external moments (see §2.2) and have been expressed 
as a part of the orbital period (AP). Errors were estimated according to the accuracy of the data and the quality of the light 
curve covering. 



shallow (Am vis w 0.8) and probably the shortest one - 
its total duration was about 22 days. The last observed 
1997 (E = 8) eclipse was, in contrary deeper than usual. 
However, its shape and the total duration were similar to 
those of 1975 and 1980. 

2.2. The characteristic phases of the minimum 

In general we can note that 1) the shape of the eclipse 
may change from the deep and narrow minimum (like in 
1964) to a shallow and very wide minimum (1969, 1986), 
2) a flat-bottom phase may occur which is characteristic 
for annular eclipses, 3) wide wings are often visible and 4) 
there is a persistent asymmetry of the light curve between 
the descending and the ascending branches of the minima. 

Although the light curve of EE Cep shows these 
changes, it is possible to single out the most typical shape 
of the minimum. MG suggested that two eclipses from 
1975 and 1980 are most typical. Basing on the data from 
these two eclipses they construct the representative (" tem- 
plate") eclipse shown in Fig. |3J There are some typi- 
cal features of the minimum which are also observed in 
other eclipses of EE Cep. The most important is the 
presence of six characteristic moments during an eclipse. 
Moments 1st, 2nd, 3rd and 4th correspond to contacts 
observed in normal eclipsing binaries, while moments l~st 
and 4 + th denote the start and the end of the phase of 
semi-atmospheric eclipse, respectively. Another feature is 
the asymmetry of the light curve between the ascending 
and descending branches of the minimum. Such asymme- 
try is clearly visible in eclipses of epochs 4, 5 and 8 and 
may be supposed to occur in the eclipse of epoch 6. The 
asymmetry seems to be caused by the presence of the 
sloping-bottom phase (between 2nd and 3rd contacts - 
Fig.EJ) which influences the shape of the ascending branch. 
This picture is supported by the fact that during eclipse 
of epoch 3 four distinct contacts were observed and the 
bottom was flat. Thus the part of the light curve between 
the 2nd and 3rd moments is not part of the egress. 

The characteristic shape and phases of the eclipse are 
independent of the photometric band. Over a broad spec- 



lst 1st 2nd 3rd 4th 4th 




10 20 30 40 50 60 

days 



Fig. 3. The characteristic shape of the EE Cep's minima 
(E=4 and E=5 joined, se for details MG). In general six 
moments can be distinguished in most of minima. The 
wings corresponding to semi-atmospheric eclipse are ap- 
parent between Pst-lst and 4-4 + th moments. The "true" 
eclipse can be defined as the part of the light curve be- 
tween the 1st and 4th moments, with a sloping-bottom 
between 2nd and 3rd contacts. 

tral range, from U until /, the eclipse seems to be almost 
grey (MG). During three eclipses (E = 4, 5, 6) the depth 
in the B band was on average 0.07 mag deeper than in 
the V band, possibly due to a selective extinction of the 
eclipsing body matter. 

2.3. Depth - duration relation 

Table ^ presents the total duration and the depth of 
EE Cep's eclipses estimated from data in Figs. ^ an d 121 
Figure 21 presents the relation between the duration and 
the depth of eclipses in B/m pg . For epochs 6 and 7 we 
have only V/vis estimates so we applied a correction of 
Am = +0.07 mag to the depth of these two eclipses. 
It is worth noting that there is a quite distinct correla- 
tion between these two parameters - longer eclipses cor- 
respond to shallower minima. Also, from the inspection 
of Figs H an d [U one can note that the shallowest minima 
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Fig. 4. The depth-duration relation for EE Cep's eclipses. 
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Fig. 5. Time dependence relation for the depth of the 
eclipses. 

have the most asymmetric wings. The minimum of Epoch 
7 seems to be an exception to the rule. However, the es- 
timates of AP and AV for this eclipse are based mainly 
on inconsistent visual observations, this may reflect the 
presence of a more complex relation than simple linear 
dependence. 

The value of AB/m pg depths versus the number of 
epoch is presented in Fig. Its evolution is rather com- 
plex, but a possible period of nine epochs (~ 50 years) can 
be suggested. Because it is almost equal to the timespan 
of all the historical data, new observations are needed to 
confirm its existence. The duration of the eclipses, which 
is a more difficult parameter to measure, do not show such 
clear time evolution. 

3. The model of disc eclipses 

The model presented by MG assumes that the invisible 
component is an optically thick, dark disc inclined to the 
orbital plane - a "cousin" of the disc in e Aur system 
(e.g. Carroll et al. 1991). The spin axis precession of the 
disc changes two angles: 1) inclination of the disc to the 
line of sight (which causes changes of the projected area 
of the disc) and 2) inclination of the projected disc to the 



orbital movement direction (tilt). As a result we observe 
varying depth of the minima as it is shown in Fig. 4 in 
MG. Another kind of a precession - the nodal regression 
- was recently proposed also for two relatively close eclips- 
ing binaries: SS Lac (Torres 2001) and V907 Sco (Lacy et 
al. 1999) as a explanation of a varying eclipse depth. In 
both cases there are precession motion (timescale of hun- 
dred years) of the whole short-period (several days) binary 
orbit caused by an interaction with a third body moving 
on a wider orbit with a period of hundreds days. In the 
case of the long period binary EE Cep, however, we ob- 
serve only the precession of the disc itself. 

During our quantitative analysis and after doing some 
simple numerical tests, we have found that the most prob- 
able reason for the observed asymmetries was not the 
asymmetry of the disc itself but the tilt of the projected 
disc to the direction of motion. The tilt produces typical 
asymmetry of the minima as the result of the sloping- 
bottom between inner contacts (Fig. . The wide wings 
between 1st and 1st, and between 4th and 4+th may be 
caused by the outer, semitransparent part of the disc. The 
"true" eclipse is caused by the opaque, inner part: the 
four characteristic moments lst-4th may be easily under- 
stood as contacts between the outer and inner parts of the 
tilted disc with the eclipsed star - see Fig. and During 
the extreme, rare, edge-on and small tilt configuration the 
eclipses should have a flat bottom phase as was observed 
in 1969. The possible period of nine epochs (~ 50 years) 
may reflect precession of the disc. 

4. Spectroscopic behaviour 

There are very few spectroscopic observations of EE Cep 
described in the literature. Herbig l|1960[) gave its spec- 
tral type as B5:ne/3. Bruckner |1976) noted that the 
emission profile in H a is not typical of Be stars and 
during the eclipse of epoch 4 there was a considerable 
strengthening of this emission in respect to the contin- 
uum. However, his observations were done in very low 
resolution (200 A/mm). Baldinelli et al. I|1981|l did not 
record any changes in the spectrum during the following 
eclipse (E=5), in spite of the slightly higher resolution (60 
A/mm). Moreover, they noted broad absorption in the 
Balmer series from H a to H% with strong, broad (H a ) 
and narrow {Hp) emission lines superimposed. This spec- 
trum was estimated as a typical Be star of B5 III spectral 
type, which agrees with the B — V, U ~ B colours analysis 
(MG). Thus, the connection between Balmer emission and 
eclipsed (Be star) or eclipsing (disc) body is not clear. 

We have secured several spectra covering the H a and 
the sodium doublet Na I regions of EE Cep. The spectra 
were taken at NAO Rozen in Bulgaria on a 2-m telescope 
equipped with a coude-spectrograph + LN2-cooled CCD 
camera at resolution 0.4 A. Unexpected large changes of 
the H a profile can be noted outside the eclipse phases on a 
time scale of several months to years (Fig. 01. The H a pro- 
files are strong and broad (~ 600 km/s), show one or two 
peaked emission with variable asymmetries. Two profiles 
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Fig. 6. The schematic explanation of the characteristic 
moments observed in the light curve of EE Cep during 
minima via the disc model. An 'edge-on' position for the 
disc was assumed for clarity. The B5 primary is marked 
by a large circle and semitransparent parts of the disc by 
a light shadow. Because of the disc tilt during the egress 
the wing may be much shorter or even unobserved - like 
during eclipse of epoch 8. 

obtained during two consecutive days are almost identical, 
however a red-shifted bulge around 60 km/s seems to be 
stronger in the later spectrum. The sodium doublet D!D 2 
shows strong single-component, probably interstellar ab- 
sorptions. However, its heliocentric radial velocity proba- 
bly changed from V r — —36 km/s in September 1998 to 
V r = — 24 km/s in November 2002. Also the equivalent 
width of D1D2 doublet seems to be variable, decreasing 
about 10% between September 1998 and November 2002. 
The mean value of the Di equivalent width EW = 0.57 A 
corresponds to Eb-v — 0.46 using the calibration given by 
Munari & Zwitter l|1997[l . which is fairly consistent with 
the reddening derived from UBV photometry (MG). 

5. The next eclipse: 2003 (E = 9) 

The moment of the mid-eclipse according to the elements 
given in MG is JD 2452794 (3 June 2003). As shown in 
Fig. 2| the eclipse can start even about 30 days before 
mid-eclipse. The reason for it is the presence of the wide 
wing during ingress which can be usually observed, thus 
the photometric observations should start at least 40 days 
before the expected mid-eclipse. The brightness of the star 
should decrease by at least 1 mag in B, V filters, but the 
observers should keep in mind that this decrease may be 



much larger: up to 2 mag! Also, a significantly shallower 
(as epochs 3 and 7) minimum can occur. If the periodicity 
of nine epochs suggested by Fig. |S] is real we can expect 
an eclipse similar to that of E=0, i.e. narrow and deep. 
Near-far infrared photometry during the eclipses will be 
very important as it would give the answer to what is the 
temperature of the companion. Also, good quality B, V or 
R light curves from the whole eclipse will be valuable to 
make numerical tests of the disc model. 

The spectroscopic observations appear to be very im- 
portant. We should expect the discovery of a shell spec- 
trum produced by the gaseous, semitransparent part of 
the disc like in e Aur (Ferluga & Mangiacapra 1991). 
Moreover, such observations would give us an additional 
clue to the answer of the origin of the variable H a line. 
The systematic observations would be very useful to ob- 
tain the spectroscopic orbit of the visible B5 component. 

The eclipsing body in EE Cep system seems to be 
unique. The confirmation (or rejection) of the disc hy- 
pothesis can be attained only by detailed photometric and 
spectroscopic study of the incoming eclipses. 
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